February 8, 200£ 



Second order additive invariants in elementary cellular 

automata 

o 

C\l ; Henryk Fuks 

D , Department of Mathematics 

Brock University 
St. Catharines, Ontario L2S 3A1, 
Canada 
hfuks@brocku. ca 



o 
u 

d 



Abstract 



We investigate second order additive invariants in elementary cellular automata 
rules. Fundamental diagrams of rules which possess additive invariants are either linear 
| or exhibit singularities similar to singularities of rules with first-order invariant. Only 

rules which have exactly one invariants exhibit singularities. At the singularity, the 
current decays to its equilibrium value as a power law t a , and the value of the exponent 
£NJ \ a obtained from numerical simulations is very close to —1/2. This is in agreements 

with values previously reported for number-conserving rules, and leads to a conjecture 
that regardless of the order of the invariant, exponent a seems to have a universal value 
of 1/2. 



1. Introduction 



> 

Cellular automata (CA) are often described as systems of cells in a regular lattice updated 
synchronously according to a local interaction rule. An interesting subclass of CA consists 
of rules possessing an additive invariant. The simplest of such invariants is the total number 
of sites in a particular state. CA with such invariant, often called "conservative CA" or 
"number-conserving CA" , generated a lot of interest in recent years PH El HH IH E] • Number- 
conserving CA can be viewed as a system of interacting and moving particles, where in 
the case of a binary rule, l's represent sites occupied by particles, and O's represent empty 
sites. The flux or current of particles in equilibrium depends only on their density, which is 
invariant. The graph of the current as a function of density characterizes many features of 
the flow, and is therefore called the fundamental diagram. 

For a majority of number- conserving CA rules, fundamental diagrams are piecewise- 
linear, usually possessing one or more "sharp corners" or singularities. There exist a strong 
evidence of universal behavior at singularities, as reported in [7] . 
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Since number- conserving CA are simplest rules with additive invariants, it would be 
interesting to consider higher order invariants and CA rules with such invariants. In 1991, 
Hattori and Takesue performed extensive study of additive invariants in discrete-time lattice 
dynamical systems jHj, not necessarily restricted to CA. They derived very general existence 
conditions, and applied them to elementary CA rules as well as reversible CA. In this paper, 
we will use their results to study second-order invariants in elementary rules, focusing mainly 
on fundamental diagrams. 

2. Number-conserving cellular automata 

In what follows, we will assume that the dynamics takes place on one-dimensional lattice of 
length L with periodic boundary conditions. Let Si(t) denote the state of the lattice site % 
at time t, where % G {0,1, . . . , L — 1}, t G N. All operations on spatial indices % are assumed 
to be modulo L. We will further assume that Sj(£) G {0, 1}, and we will say that the site % 
is occupied (empty) at time t if Sj(i) = 1 (sj(t) = 0). 

Let I and r be two integers such that I < < r, and let n = r — I + 1. The set 
{si + i(t), Si + i + i(t), . . . , Si +r (t)} will be called the neighbourhood of the site Sj(t). Let / be a 
function / : {0, l} n — > {0, 1}, also called a local function The update rule for the cellular 
automaton is given by 

Si(t + 1) = f(s i+ i(t), s i+l+1 (t), s i+r (t)). (1) 

In jS], the concept of additive invariant for CA has been introduced. Let a be a non- 
negative integer, and let £ = £(xo, 

function of a + 1 variables taking values 
in R, We say that £ is a density function of an additive conserved quantity if for every 
positive integer L and for every initial condition (so(0), si(0), . . . , sl_i(0)) G {0, 1} L we have 

L-1 L-1 

£ £(*(*), ■ ■ ■ , *+«(*)) = £ + 1), + !),•••, + 1)) (2) 

i=0 i=0 

for all ( G N, For simplicity, if the above condition is satisfied, we will say that £ is an 
additive invariant of /. It is often more convenient to write (J2J) using the function G defined 

as 

G(yXQ, X\, . . . , 3^a+n— l) £(/(-^0) -^li • • • i •^n—l)i f\%lj -^2i ■ ■ ■ > ^n); • • • j fy^at -^a+li • ■ ■ 1 -Ka+n—l))- 

(3) 

With this notation, £ is an additive invariant of / if 

L-1 L-1 

^ ^ G{Xi, . . . , Xj_|_ Q ,_|_ n _x) ^ ^ •^i+l) • • • ) -Ei+a) (4) 

i=0 i=0 

for every positive integer L and for all x , X\, . . . , £l-i G {0, 1}. 

In recent years, many authors studied the case of the simplest additive invariant, with 
a = and £(xq) = xq. For this invariant, the equation (J2J) becomes 

L-1 L-1 

£ Si (0=5>(* + l), (5) 

i=0 i=0 
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which means that the CA rule posessing this invariant conserves the number of sites in state 1. 
Such rules are often referred to as number- conserving rules. Among elementary CA, i.e. those 
with I — — 1 and r = 1, there are only five number-conserving rules. Three of these are trivial, 
namely the identity rule 204 and two shifts 170 and 240. Two remaining rules, 184 and 226, 
are equivalent under the spatial reflection. Rule 184, which is a discrete version of the totally 
asymmetric exclusion process, has been extensively studied [TUl HH H2 Ed 113 113 EE] , and 
many rigorous result regarding its dynamics have been established. 

Hattori and Takesue [Sj established a very general result which we will write here in a 
somewhat simplified form, taking into account that this paper is concerned with binary rules 
only. 

Theorem 1 (Hattori &; Takesue '91) Let £,(x ,Xi, • • • ,x a ) be a function of a + 1 vari- 
ables. Then £ is a density function of an additive conserved quantity under the time evolution 
of cellular automaton rule (QJ) if and only if the condition 

G^Xq, Xi, . . . , X a+n _i) — £(x_i, . . . , X a ^{) = J (Xq, Xi, . . . , X a+n -2) ~ J( x li %2i ■ ■ ■ > x a+n-l) 

(6) 

holds for all Xq, X\, x a+n -i G {0,1}, where the quantity J, to be referred to as the 
current, is defined by 

J(x ,x 1 ,. . .,x a+n - 2 ) = - 22 Gr ( ' ' • • • >°> x o,£i, • • • ,Xi)+ 2j ^(0, 0, . . . ,0,x ,x 1 , . . . ,Xi-t) 

i=0 i=-l-n+2 

(7) 

The following convention is used in the definition of J: 

a+l—i 

£('0, 0, . , 0, x , xi,..., Xi-i) = ^(xi-a-i, Xi- a , . . . , Xi-i) ifi>a+l, (8) 

and 

a+l—i 

e(5J A 3,a; o ,a; 1) ... 1 a: i _ 1 ) = e(0 l 0,...,0) if i < 0. (9) 

The equation © can be interpreted in a similar way as a conservation law in a continuous, 
one dimensional physical system. In such system, let p(x, t) denote the density of some 
material at point x and time t, and let j(x, t) be the current (flux) of this material at point x 
and time t. A conservation law states that the rate of change of the total amount of material 
contained in a fixed domain is equal to the flux of that material across the surface of the 
domain. The differential form of this condition can be written as 

I -I- « 10 > 

Since in our case £ is the the density of an additive conserved quantity, the left hand side of 
© is simply the change of density in a single time step, so that © is an obvious discrete 
analog of the current conservation law (|10j) with J playing the role of the current. 
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Let us now assume that the initial configuration has been generated from some translation- 
invariant distribution p. We define the expected value of £ at site i as 

p(i, t) = [Z( Si (t), s i+1 (t), . . . , s i+a (t))} . (11) 

Since the initial distribution is i-independent, we expect that p(i,t) also does not depend 
on i, and we will therefore define p(t) = p(i,t). Furthermore, since £ is density function of 
a conserved quantity, p(t) is t- independent, so we define p = p(t). The expected value of 
the current J(s i+ i(t), . . . , s i+r _i(t)) will also be ^-independent, so we can define the 

expected current as 

j(p,t) = E fl (j(s i+ i(t),s i+ i +1 (t),...,s i+r - 1 (t))). (12) 

The graph of the equilibrium current j(p, oo) = lim^ooj(p, t) versus the density p is 
known as the fundamental diagram. 

3. Number-conserving nearest-neighbour rules 

In order to illustrate the theorem of the previous section, we will first consider the case of 
number-conserving nearest-neighbour rules, i.e., a = and (,(x ) — x , I — —1, n = 3. 
Condition © becomes 

G(x , xi, x 2 ) - i{xi) = J (x , xi) - J (xi,x 2 ), (13) 

where 

J(x , Xl ) == -G(0, 0, x ) - G(0, x , xi) + £(0) + £(x ) (14) 

Obviously, G(x ,xi,x 2 ) = f(x ,xx,x 2 ), thus the current becomes J(x ,Xi) = — f(0, 0, Xq) — 
f(0,x ,Xi) — xq, and the conservation condition takes the form 

f(x ,x 1 ,x 2 ) -xi = J(x ,xi) - J(x 1 ,x 2 ). (15) 

As mentioned earlier, rule 184 and its spatial reflection are the only non-trivial elementary 
CA rules satisfying (fT3j). For rule 184, / is defined by f(x , xi, x 2 ) = x — x Xi + X\X 2 and 
the current can be written as J(x ,Xx) — x (l — Xi). It is possible to show [13] that the 
equilibrium current for this rule is given by 

■( _ / P if P < 1/2, , lg s 
J\.P,°°) I 1 _ p otherwise. ^ ' 

Since number-conserving CA rules conserve the number of occupied sites, we can label 
each occupied site (or "particle") with an integer k e Z, such that the closest particle to 
the right of particle k is labeled k + 1. If Ukif) denotes the position of particle k at time t, 
the configuration of the particle system at time t is described by the increasing bisequence 
{yk(t) }j^L-<xr We can then specify how the position of the particle at the time step t + 1 
depends on positions of the particle and its neighbours at the time step t. For example, for 
rule 184 one obtains 

y k (t + 1) = y n (t) + mm{y k+1 (t) - y k (t) - 1, 1}. (17) 
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Equation (JT7j) is sometimes referred to as the motion representation. The motion representa- 
tion is analogous to Lagrange representation of the fluid flow, in which we observe individual 
particles and follow their trajectories |17j . It turns out that the motion representation can 
be constructed for arbitrary number-conserving CA rule by employing algorithm described 
in [IE!. 

4. Second-order invariants 

In what follows, we will referr to the number of variables of £ as the order of the invariant, 
equal to a + 1. Since the invariant of a = and corresponding fundamental diagrams have 
been extensively studied, we will explore the case of a = 1, i.e., second order invariants, 
using the method of [8 . 

The arguments Xq, X\ of the density function take values in the set {0, 1}, and therefore 
£ can be defined in terms of four parameters 

£(0, 0) = coo, £(0, 1) = cox, f (1, 0) = do, £(1, 1) = di, (18) 

where c o, c i, do> di £ This can be also expressed as 

£(x , xi) = c 00 (l - x )(l - £i) + c i(l - x )x 1 + c w x (l - xx) + cuXqXx 
— c oo + ( c io — Coo)a^o + ( c oi — Coo)a^i + (ox) — On — do + dij^o^i- 

The constant term does not bring anything new, so we can set Coo = 0. Moreover, note that 
for any function g(x) and any xq, Xx, ■ ■ ■ , Xl-x £ {0, 1} we have 

L-1 L-1 

/] £(Xj, Xi + x 

i=0 i=0 

which means that if £(xq,Xi) is a density function of some conserved additive quantity, then 
F(xq, Xi) = £(x , Xx)+g(xo)—g(xx) is also a density function of a conserved additive quantity. 
To remove this ambiguity, we will require that £(0, Xi) = 0, similarly as done in This 
yields c i = 0, and we are left with £ depending on two parameters only 

£(x , xx) = c 10 xo + (di - c 10 )x xx. (20) 

Defining a\ = —do, a 2 — di — do, we arrive at the final parameterization of £ 

£(x ,a;i) = axx + a 2 x xx, a 1 ,a 2 e'R. (21) 

For a = 1, I = —1, and n = 3, eq. (jOJ becomes 

G(xo,Xi,x 2 ,x 3 ) - £(xi,x 2 ) = J(x ,xi,x 2 ) - J(xi,x 2 ,x 3 ), (22) 

where 

G(a; ,a;i,x 2 ,a;3) = i(f(xo,x 1 ,x 2 ),f(x 1 ,x 2 ,x 3 )) 

J{xq,xx,x 2 ) = -G(0,0,0,x ) - G(0,0,x ,xi) - G(0,x o ,xx,x 2 ) 

+ e(o,o) + £(o,x ) + e(x ,x 1 ). 
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Rule number 


£(z ,zi,2r 2 ) 


J(x ,xi) 


12 


Xo — XqX± 


— XoXi 


14 


Xq - XqXi 


— XqX\ 


15 


Xq — XqX\ 


-XqX 1 


34 


Xq — XqX\ 


-Xi 


35 


xq - x x 1 


-Xi 


42 


Xq — XqX± 


—X\ + X\X 2 — XqX\X 2 


43 


Xq — XqX± 


—X\ + X\X 2 — XqX\X 2 


51 


Xq — XqX± 


-Xi 


140 


Xq — XqX± 


-XqXi + XqXiX 2 


142 


Xq — XqX± 


-XqXi + XqXiX 2 


200 


XqXi 






Table 1: Density of the invariant £ and the current J for all non-trivial elementary CA with 
second order invariants. 

Since £(0, 0) = £(0, x ) = 0, the formula for current simplifies to 

J (x , x 1 ,x 2 ) = -G(0, 0, 0, x ) - G(0, 0, x , xx) - G(0, x , x lt x 2 ) + £(ar , x{). (23) 

For a given elementary CA rule f(xQ,Xi,x 2 )i one can write eq. (J22|) for all 2 4 combinations 
of values of the variables xq,xi,x 2 ,x^ G {0,1}, thus obtaining an overdetermined linear 
system of 16 equations with two unknowns ai,a 2 . This system is homogeneous, therefore 
the solution, if it exists, is not unique. That is, if (ai, a 2 ) is a solution, then (cai, ca 2 ) is also 
a solution for any c G EL We will normalize the solution so that the first non-zero number 
in the pair (ai, a 2 ) is set to be equal to 1. 

Solving these equations for all "minimal" CA rules 1 , one finds that for most CA rules 
solutions do not exist. Remaining CA rules can be divided into two classes. The first class 
contains rules 204 and 170, and for these rules, any pair (ai, a 2 ) is a solution. We will not be 
concerned with these rules, since they exhibit trivial dynamics. The second class consists of 
10 rules for which a unique solution exists (up to the normalization described earlier). These 
rules are 12, 14, 15, 34, 35, 42, 43, 51, 140, 142, and 200, as reported in 0. 

Table 1 shows the density function £ and the current J for all of them. The formulas 
for the current have been obtained using the HCELL C++ library for cellular automata 
developed by the author. 

5. Fundamental diagrams 

In order to construct fundamental diagrams for rules of Table 1, we first note that the current 
for rule 200 is identically equal to zero, thus the equilibrium current j(p, oo) = 0. The graph 
of j(p, oo) vs. p for this rule is, therefore, not interesting. 

1 Elementary CA rules fall into 88 equivalence classes with respect to the group of transformations gen- 
erated by the spatial reflection and the Boolean conjugacy. Minimally-numbered element of each class are 
known as "minimal rules" . 
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For all other rules, the density of the invariant is given by the same function £(xq,xi) = 
Xq — XqX\. This means that rules 12, 14, 15, 34, 35, 42, 43, 51, 140, and 142 conserve 
the number of blocks "10" in the configuration. In order to construct their fundamental 
diagrams, we have to be able to create an initial configuration with a given number of pairs 
"10" . Construction of a configuration of length L with exactly m pairs "10" can proceed 
according to the following algorithm. We start with an array of L integers, {sj(O)}^, 1 . 

1. Set Si(0) = for all % = 0, 1, . . . , L - 1. 

2. Place the symbol "C" at randomly selected site of the array. Then place another symbol 
"C" at another site randomly selected among all remaining empty sites. Repeat this 
procedure until you place exactly 2m symbols "C" . 

3. Let x = 0. Starting from % — 0, traverse the array filling it with x values. Every time 
when you encounter C, set x := 1 — x. Stop when you reach the end of the array. 

The average density of the invariant £(xq,Xi) = x — x Xi for the configuration obtained 
with the above algorithm will be 

1 i_1 m 
P™ = 1^^ i (t),s i+1 (t)) = j, (24) 

i=0 

and it will be independent of t. We can also define average current at time t for a configuration 
with average density p av of the invariant as 

1 i_1 

jav(Pav,t) = -22j(Si(t),Si +1 (t),Si +2 (t)). (25) 
i=0 

Graph of j a v{Pav,t) vs. p av for very large t will approximate the graph of j(p, oo) (given by 
eq. H2J) vs. p, i.e., the fundamental diagram. 

For six rules from Table 1, the fundamental diagram is strictly linear, and the following 
expressions for current can be conjectured based on numerical experiments. 



Rule 12 


j(p,oo) 


= 




Rule 15 


j(p,oo) 


= p 


-1/2 


Rule 34 


j(p,oo) 




P 


Rule 42 


j(p,oo) 




P 


Rule 51 


j(p,°°) 




1/2 


Rule 140 


j(p,oo) 


= 





The remaining four rules are more interesting, as they exhibit singularities in fundamental 
diagrams, as shown in Figure 1. It is remarkable that singularities are present only in 
fundamental diagram of those particular rules. In [Sj, authors searched for invariants of up 
to seventh order for all "minimal" elementary CA. According to the table published in their 
paper, rules 14,35,43, 142 do not have any other invariant except £(xo,Xi) = xq — xqXi, in 
contrast to remaining rules of Table 1, which also posses other higher order invariants. It 
seems that singularities in the fundamental diagram can appear only in rules which have 
only one invariant, just like rule 184, which only has first order invariant £(xq) = Xq. 
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Figure 1: Fundamental diagrams for rules 14, 35, 43, 142. Singularities are denoted by green 
circles. Diagrams have been obtained using 50000 lattice sites after 50000 iterations. 



S 







0.1 



0.15 



0.2 



0.25 



0.3 



0.35 



0.4 



P 



Figure 2: Decay time vs. density of the invariant for rule 43. 



6. Convergence to equilibrium 

In number-conserving cellular automata, singularities of the fundamental diagram exhibit 
critical behavior. This can be illustrated by introducing the the decay time defined as 



If the decay of j(p,t) toward its equilibrium value j(p, oo) is of power-law type, the above 
sum diverges. For all rules in Figure 1, we have performed computer simulations to estimate 
r. The value of r has been estimated by measuring j a v(p,t) for t = 0, 1, . . . , 1000, and 
truncating the sum (}2T)|) at t = 1000. Figure 2 shows a typical graph of r as a function of 
p, obtained for rule 42. Comparing Figures lc and Figure 2 we clearly see that r diverges 
at the critical point of rule 42, which occurs at p = 0.25. We will denote this value by p c . 
Assuming that |j(p c , oo) — j(p c ,t)\ ~ t~ a , we have determined the exponent a as the slope 
of the straight line which best fits the logarithmic plot of |j(p c , oo) — j(p c ,t)| vs. time t. 
Example of such a plot, again for rule 42, is shown in Figure 3. Table 2 shows values of the 
exponent a for critical points of all four rules of Figure 1. 

Exponent a is known to be equal to exactly 1/2 for rule 184 and its generalizations, and 
rigorous proof of this fact exists [13j. Extensive numerical experiments support the conjecture 
that for all rules with first-order invariant the value r = 1/2 is universal, in the case of both 
piecewise linear [H] and nonlinear jB] fundamental diagrams. Table 2 provides evidence that 
a more general conjecture may be valid: regardless of the order of the invariant, exponent a 




(26) 



t=o 
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Figure 3: Logarithmic plot of \j(p c ,oa) — j(p c ,t)\ as a function of time for rule 43. Data 
points (+) represent computer simulations, while the dashed line represents the best fit. 



Rule number 


(pc,j(pc,oo)) 


a 


14 


(1/4, -1/4) 


-0.504 


35 


(1/3,-1/3) 


-0.472 


43 


(1/4, -1/4) 


-0.492 


142 


(1/4, -1/4) 


-0.502 



Table 2: Values of the exponent a at the critical point for elementary CA rules with second- 
order additive invariant. 
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Figure 4: Spatiotemporal pattern generated by rule 184. Initial configuration is represented 
by the top row, where black squares represent l's, and white spaces represent O's. Consecutive 
configurations are plotted as consecutive rows. Two types of defects of "A" and "B" type 
are visible, annihilating in the circled spot. 

seems to have universal value of 1/2. In the next section we will offer some justification for 
this conjecture for rules with second-order invariants. 

7. Dynamics of localized structures 

In rule 184, the power-law convergence of the current toward its equilibrium value is related 
to the dynamics of this rule, which resembles ballistic annihilation. The spatiotemporal 
patter generated by rule 184 can be understood as propagation of two types of localized 
structures, shown in Figure 4. These two types of structures, marked with letters "A" and 
"B", propagate in opposite directions and annihilate upon collision. At the critical point, the 
number of "A" defects in the initial configuration is the same as the number of "B" defects, 
and it takes long time for all of them to disappear, hence the "critical slowing down", or 
power-law convergence is observed. Detailed analysis of this process [T3] leads to the exact 
formula for the current j(p,t), which in the limit of large t and using de Moivre-Laplace 
limit theorem leads to j(p c , oo) — j(p ci t) ~ t -1 / 2 . Here, by f(t) ~ g(t) we mean that 
lim^oo f(t)/g(t) exists and is different from 0. 

Dynamics of rules 14, 35, 43, 142 resembles rule 184 very strongly, as can be seen in 
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Figure 5, which shows spatiotemporal patterns at the critical point for all four rules. In all 
four cases, localized propagating structures moving in opposite directions and annihilating 
upon collision are visible. In fact, for two of these rules, it is possible to establish direct 
relationship with rule 184. In order to do this, we will define superposition of two rules as 



If ho f — ho g 7 then following 19 we will say that g is a transform of rule / by h. If by fk 
we denote local function of rule k, one can show ^Hl that 



This means that there exists a local mapping (rule 60) which transforms rule 43 into rule 
184, and rule 142 into rule 226 (recall that rule 226 is the image of rule 184 under spatial 
reflection). Similarity of dynamics of rules 43, 142 to rule 184 is, therefore, not surprising. 

8. Conclusion 

We investigated second order additive invariants in elementary cellular automata rules. We 
found that fundamental diagrams of rules which possess additive invariant are either linear 
or exhibit singularities similar to singularities of rules with first-order invariant. Singularities 
can appear only in rules with exactly one invariant. At the critical density of the invariant, 
the current decays to its equilibrium value as a power law t a , and the value of the exponent 
a obtained from numerical simulations is very close to —1/2. This indicates that regardless 
of the order of the invariant, the dynamics of CA rules with invariants is very similar. 

Since rules 43 and 142 can be transformed into rules 184 and 226 by a surjective local 
transformation, it should be possible to obtain for them rigorous formulas for the expected 
value of the current at arbitrary time, similarly as it has been done for rule 184 and its 
generalizations [13j. Such formula could then used to compute the exact value of the exponent 
a. For rules 14 and 35 no such local transformation exists, nevertheless they exhibit localized 
propagating structures strikingly similar to structures of rule 184, so exact calculation of the 
current might be possible too. This problem is currently under investigation and will be 
reported elsewhere. 

Acknowledgements: The author acknowledges financial support from NSERC (Natural 
Sciences and Engineering Research Council of Canada) in the form of the Discovery Grant. 

References 

[1] M. Pivato, "Conservation laws in cellular automata," Nonlinearity 15 (2002) 
1781-1793, arXiv: math .DS/01 11014] 

[2] A. Moreira, "Universality and decidability of number- conserving cellular automata," 
Theor. Comput. Sci. 292 (2003) 711-721. 



(/ ° fiOOo, x 1 ,x 2 , x 3 , x 4 ) = f(g(x , xi,x 2 ),g(x 1 , x 2 , x 3 ),g(x 2 , x 3 , x 4 )). 



(27) 



/60 ° fi3 — /l84 ° foo, 
iW) ° /l42 = ^226 ° ./60- 



(28) 
(29) 



12 




Figure 5: Spatiotemporal pattern generated by rules 14, 35, 43, 142 at their critical points. 
Density of the invariant £ equals 1/4 for rules 14, 43, 142, and 1/3 for rule 35. 



13 



B. Durand, E. Formenti, and Z. Roka, "Number-conserving cellular automata I: 
decidability," Theoretical Computer Science 299 (2003) 523-535. 

E. Formenti and A. Grange, "Number conserving cellular automata II: dynamics," 
Theoretical Computer Science 304 (2003) 269-290. 

K. Morita and K. Imai, "Number-conserving reversible cellular automata and their 
computation-universality," Theoretical Informatics and Applications 35 (2001) 
239-258. 

H. Fuks and N. Boccara, "Convergence to equilibrium in a class of interacting particle 
systems," Phys. Rev. E 64 (2001) 016117, arX iv : nlin . CG/01010371 

H. Fuks, "Critical behaviour of number-conserving cellular automata with nonlinear 
fundamental diagrams," J. Stat. Mech.: Theor. Exp. (2004). art. no. P07005. 

T. Hattori and S. Takesue, "Additive conserved quantities in discrete-time lattice 
dynamical systems," Physica D 49 (1991) 295-322. 

J. Krug and H. Spohn, "Universality classes for deterministic surface growth," Phys. 
Rev. A 38 (1988) 4271-4283. 

T. Nagatani, "Creation and annihilation of traffic jams in a stochastic assymetric 
exclusion model with open boundaries: a computer simulation," J. Phys. A: Math. 
Gen. 28 (1999) 7079-7088. 

K. Nagel, "Particle hopping models and traffic flow theory," Phys. Rev. E 53 (1996) 
4655-4672, larXiv : cond-mat/9509075] 

V. Belitsky and P. A. Ferrari, "Invariant measures and convergence for cellular 
automaton 184 and related processes," math. PR/9811103, Preprint. 



H. Fuks, "Exact results for deterministic cellular automata traffic models," Phys. Rev. 
E60 (1999) 197-202, [arXiv: comp-gas/9902001[ 

K. Nishinari and D. Takahashi, "Analytical properties of ultradiscrete burgers equation 
and rule-184 cellular automaton," J. Phys. A- Math. Gen. 31 (1998) 5439-5450. 

V. Belitsky, J. Krug, E. J. Neves, and G. M. Schutz, "A cellular automaton model for 
two-lane traffic," J. Stat. Phys. 103 (2001) 945-971. 

M. Blank, "Ergodic properties of a simple deterministic traffic flow model," J. Stat. 
Phys. Ill (2003) 903-930. 

J. Matsukidaira and K. Nishinari, "Euler-Lagrange correspondence of cellular 
automaton for traffic-flow models," Phys. Rev. Lett. 90 (2003) art. no.-088701. 



14 



[18] H. Fuks, "A class of cellular automata equivalent to deterministic particle systems," in 
Hydrodynamic Limits and Related Topics, A. T. L. S. Feng and R. S. Varadhan, eds., 
Fields Institute Communications Series. AMS, Providence, RI, 2000. 
larXiv : nlin . CG/02070471 

[19] N. Boccara, "Transformations of one- dimensional cellular automaton rules by 
translation-invariant local surjective mappings," Physica D 68 (1992) 416-426. 



15 



